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Abstract 

The probability distribution of log-returns of financial time series, sampled at high fre- 
quency, is the basis for any further developments in quantitative finance. In this letter, we 
present experimental results based on a large set of time series on futures. Then, we show 
that the t-distribution with v ~ 3 gives a nice description of almost all data series. This 
appears to be a quite general result that stays robust on a large set of any financial data as 
well as on a wide range of sampling frequency of these data, below one hour. 



1 Introduction 



Returns in financial time series are the most fundamental inputs to quantitative finance. To a 
certain extend, they provide some insights in the the dynamical content of the market. In this 
letter, we consider several financial series on futures, using always a five minutes sampling. 
Each future contract is characterized by a price series S(t), from which we extract the log- 
returns x(t) = log[S(t + 1)/S(t)}. The analysis is then driven on these log-returns x(t). 

From standard quantitative analysis IBBl, we know that the distribution of log-returns x(t), 
namely P(x), can be written quite generally as 

P{x) = ^exp(~w(x)/2) (1) 

where w(x) is an objective function and Z a normalization factor. In particular, it can be shown 
easily that w(x) can be derived by minimizing a generating functional F[w(x)], subject to some 
constraints on the mean value of the objective function. It reads 

F = /<fcP(*)[.o g P(*) + M*)/Z>-A] (2) 

where A is an arbitrary constant. 

In addition, the expression given in Eq. (OQ) for the probability distribution can also be seen 
as the outcome of an equation of motion for x(t). From Eq. (OQ) and ©, we can express the 
stochastic process x(t) as a Markovian process of the form 03-IH 

(IT 

- = f( x ) +g ( x )e(t) (3) 

where e(t) is a Gaussian process satisfying < e(t)e(t') >= D5(t — t') and < e(t) >= 0. In Eq. 
©, functions / and g depends only on x{t). Adopting the Ito convention iBBl, the distribution 
function P(x,t), associated with this equation of motion (Eq. ©), is given by the following 
Fokker-Planck equation 

d - E ¥ 1 = |^ 2 (*)p(m)] - ^[mp(x,t)} ( 4) 

From Eq. ©, we can finally extract the stationary solution for P(x) in the form of Eq. (QQ) 



P(x) = — exp 

Whether Eq. ©, © and © can be related to real data on financial markets is not granted. 
Therefore, we need to compare predictions derived from these equations to real data. As men- 
tioned above, we use financial time series on different futures, using a five minutes sampling. 
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2 Parameterizations of P(x) 



In Eq. ©, © and ©, functions / and g are not specified and any choice can be considered. 
Obviously, only specific choices will have a chance to get a reasonable agreement with real 
data. For example, let us consider three cases: 

(i) If f(x) = —x and g(x) = 1, we obtain P(x)Z = exp(—x 2 /D), and thus we predict a 
Gaussian shape for the log-returns distribution. 

(ii) In the more general case where f(x) = Ag^| and g is not constant, we obtain 

P{x)Z 



g 2(l-X/D) 



and thus we predict non-Gaussian shape for the log-returns distribution. 

(iii) Let us specify the case (ii). Introducing a constant v and defining the two functions / and 

g as f(x) = vO- + 7r) and = (1 + v)> we S et 

P(x)Z = 1 . (6) 

V > (l + ^)(v+l)/2 V 

In this scenario P(x) follows the so-called t-distribution. It depends on one parameter v 
to be fitted on real data, for normalized log-returns. Let us notice that Eq. © is equivalent 
to the q-exponential form of Ref . j51[6l 

where we have conserved the notations of Ref. J51[6l. Obviously, Eq. © and © are 
directly related by {y + l)/2 = l/(g — 1). In particular, v = 3 is equivalent to g = 1.5. 



3 Experimental Analysis of P(x) 

In Fig. [H we present the log-returns x (five minutes sampling) for a large set of futures. To make 
the comparison, we have scaled x for all futures to the volatility of the DAX future (FDAX). 
On the left hand side (LHS) of Fig. [T] we observe that futures on DAX, Bund, Yen, Euro, 
Gold present the same probability distribution for x, hence P(x) is universal for all these data 
series once the volatility is normalized to the same value. On the right hand side of Fig. [TJ we 
provide comparisons with futures on commodities. We observe that data series on CL (Crude 
Oil) follows the same P(x) as FDAX, but other data series on Wheat and NG (Natual Gaz) 
exhibit some larger tails. 

We can use results developed in section |2] in order to compare with experimental distribu- 
tions P(x) of Fig. CD We use predictions exposed in cases (i) and (iii), respectively the Gaussian 
and the q-exponential forms (or t-distribution). Results are shown in Fig. [2l For data, we only 
display P(x) for FDAX. 
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In Fig. [21 we show that the q-exponential probability distribution of Eq. ©, with q = 1.5, 
gives a good description of the data. Similarly, it corresponds to v = 3 in the form of the t- 
distribution of Eq. ©. Let us notice that futures with larger tails discussed above correspond to 
smaller values of v. Also, we observe in Fig. [2]that the Gaussian approximation fails to describe 
properly the data. On the right hand side of Fig. [2l we observe that when v is increased above 
3 in Eq. ©, then P{x) stands in the middle of the correct probability density and the Gaussian 
approximation. When v tends to infinity, the t-distribution recovers the Gaussian limit. 



4 Conclusion 



From the theory point of view, the probability distribution of log-returns of financial time series, 
sampled at high frequency, can be expressed quite generally as P(x) = exp(—w(x)/D). In 
this expression, the function w(x) can be derived using the formalism of stochastic differen- 
tial equations applied to quantitative finance. Several cases have been considered in this letter. 
Then, we have shown that real data (with normalized log-returns) are compatible with a distri- 
bution of the form 

/ x 1 

P(x) oc 5— , V ~ 3. 

V ' (1 _|_ 2i)(«H-l)/2 ' 

We have exemplified this property on a sample of data sets. However, this is a quite general 
result that stays robust on a larger set of financial data as well as on a wide range of sampling 
frequency of these data, below one hour. Interestingly also, this probability distribution corre- 
sponds to the equation of motion Eq. [3l with the expressions of / and g given in case (iii). This 
can serve as a basis for a Monte-Carlo simulation of market data and subsequent risk analysis. 
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Figure 1: Log-returns x (five minutes sampling) for a large set of futures. To make the com- 
parison, we have scaled x for all futures to the volatility of the DAX future (FDAX). Left: we 
observe that futures on DAX, Bund, Yen, Euro, Gold present the same probability distribution 
for x, hence P(x) is universal for all these data series once the volatility is normalized to the 
same value. Right: we observe that data series on CL (Crude Oil) follows the same P(x) as 
FDAX, but other data series on Wheat and NG (Natual Gaz) exhibit some larger tails. 
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Figure 2: We show that the q-exponential probability distribution of Eq. ©, with q = 1.5, gives 
a good description of the data (Right). It corresponds to v = 3 in the form of the t-distribution 
of Eq. © (Left). The Gaussian approximation fails to describe properly the data. When v is 
increased, we approach the Gaussian limit. 
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